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Summary 

Background The new Middle East respiratory syndrome coronavirus (MERS-CoV) infection shares many clinical, 
epidemiological, and virological similarities with that of severe acute respiratory syndrome (SARS)-CoV. We aimed to 
estimate virus transmissibility and the epidemic potential of MERS-CoV, and to compare the results with similar 
findings obtained for prepandemic SARS. 


Methods We retrieved data for MERS-CoV clusters from the WHO summary and subsequent reports, and published 
descriptions of cases, and took into account 55 of the 64 laboratory-confirmed cases of MERS-CoV reported as of 
June 21, 2013, excluding cases notified in the previous 2 weeks. To assess the interhuman transmissibility of MERS- 
CoV, we used Bayesian analysis to estimate the basic reproduction number (R,) and compared it to that of prepandemic 
SARS. We considered two scenarios, depending on the interpretation of the MERS-CoV cluster-size data. 


Results With our most pessimistic scenario (scenario 2), we estimated MERS-CoV R, to be 0-69 (95% CI 0-50-0-92); 
by contrast, the R, for prepandemic SARS-CoV was 0-80 (0-54—1-13). Our optimistic scenario (scenario 1) yielded a 
MERS-CoV R, of 0-60 (0-42-0-80). Because of recent implementation of effective contact tracing and isolation 
procedures, further MERS-CoV transmission data might no longer describe an entire cluster, but only secondary 
infections directly caused by the index patient. Hence, we calculated that, under scenario 2, eight or more secondary 
infections caused by the next index patient would translate into a 5% or higher chance that the revised MERS-CoV R, 
would exceed 1—ie, that MERS-CoV might have pandemic potential. 


Interpretation Our analysis suggests that MERS-CoV does not yet have pandemic potential. We recommend enhanced 
surveillance, active contact tracing, and vigorous searches for the MERS-CoV animal hosts and transmission routes to 


human beings. 


Funding Agence Nationale de la Recherche (Labex Integrative Biology of Emerging Infectious Diseases), and the 


European Community’s Seventh Framework Programme project PREDEMICS. 


Introduction 

Since September, 2012, WHO has been notified of 
64 cases of Middle East respiratory syndrome coronavirus 
(MERS-CoV) infection.’ This new disease has many 
features reminiscent of severe acute respiratory syn- 
drome (SARS)—eg, predominance of respiratory symp- 
toms, airborne transmission, high case-fatality ratio, and 
caused by a virus of the Betacoronavirus genus that is 
closely related to bat coronaviruses.** The SARS-CoV 
pandemic had been preceded by outbreaks in southeast 
China between November, 2002, and January, 2003.°’ 
During this period, the virus is thought to have been 
repeatedly introduced into human populations from its 
intermediate hosts (eg, the masked palm civet) before 
adapting to interhuman transmission.** The molecular 
events leading to the mutation that enabled the SARS- 
CoV spike protein to bind with the human ACE2 receptor 
have been particularly well described." The same 
sequence of events might be occurring with MERS-CoV, 
which has been repeatedly introduced into the human 
population for more than a year (from an unknown 
animal host) and might have human pandemic potential. 
One useful indicator of virus transmissibility is the basic 
reproduction number (R,), which represents the number 


of secondary cases per index case in a fully susceptible 
population. When R, is above 1, epidemic potential has 
been reached. We aimed to adapt a recently published 
method of estimating R,” to the MERS-CoV outbreaks 
and compare the results with similar findings obtained 
for prepandemic SARS. 


Methods 

Data sources 

We retrieved data for MERS-CoV clusters and clinical 
disease progression from various sources, including the 
WHO summary’ and subsequent reports, and published 
descriptions of cases in the UK,‘ Jordan,” France,” and 
Saudi Arabia.“ As of June 21, 2013, 64 laboratory- 
confirmed MERS-CoV cases have been reported. 

The earliest two confirmed cases appeared in a cluster 
of 13 individuals with respiratory symptoms in a hospital 
in Zarqa, Jordan, in April, 2012. Diagnosis was done 
retrospectively, by analysis of samples from these patients; 
two patients, those who died, tested positive for MERS- 
CoV.” Since then, other cases have occurred in and 
outside the Middle East (figure 1). 19 of these cases were 
regarded as sporadic, including three cases in individuals 
who were admitted to hospital in the UK and Germany. 
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Figure 1: Map of Middle East respiratory syndrome coronavirus clusters included in the analysis 
Cluster sizes are shown in bold white text. We used a blue background if their location could be established within the country of origin and a red background 
otherwise. Each arrow corresponds to travel of one patient with Middle East respiratory syndrome coronavirus infection outside the Middle East, where they caused 


secondary cases. gov=governorate. 


15 reported cases were grouped into six small clusters of 
two to four cases, including four clusters in which 
interhuman transmission occurred outside the Middle 
East (UK, France, Tunisia, and Italy). In May, 2013, a large 
cluster of 23 confirmed cases was reported in the Al-Hasa 
governorate, Saudi Arabia. An unrelated cluster of five 
cases was later reported in the neighbouring city of 
Dammam, Saudi Arabia. 

We obtained data on prepandemic SARS clusters in the 
Guangdong province of China from the scientific 
literature.” Between November, 2002, and January, 2003, 
seven clusters of one, two, three, five, seven, eight, and 
nine cases, each having one index patient, were reported. 


Data interpretation 

Our data analysis is based on the concept of transmission 
tree, defined as all cases related to one index patient. 
By contrast, the concept of cluster is based on temporal 
and geographical grouping and might include one or 
more transmission trees, depending on the number of 
index patients belonging to the cluster. 

For the branching process analysis, we assumed that 
all MERS-CoV introductions into the human population 
occurred in the Middle East. We took into account 55 of 
the 64 laboratory-confirmed cases. We excluded the nine 
cases notified in the 2 weeks preceding June 21, 2013, 
because the incubation period of MERS-CoV could be up 


to 14 days, so these cases might belong to current 
transmission clusters.* For improved modelling consis- 
tency, we also took into account seven probable cases: (1) 
one case, a patient who travelled from Saudi Arabia to 
Tunisia where two of his children, who had subsequent 
contact with him, tested positive; (2) the fourth case of 
the family cluster detected in Riyadh, Saudi Arabia, in 
October, 2012 (a patient aged 16 years who was sympto- 
matic but was not tested for MERS-CoV);* (3) cases A 
and C, reported as probable in the investigation of the 
Al-Hasa cluster;* and (4) three cases in the Jordanian 
cluster. The large number of unconfirmed cases (11 of 13) 
in the Jordanian cluster called for caution. On the basis 
of a 38% case-fatality ratio in individuals younger than 
60 years (see summary statistics in the Results section), 
and considering that the patients with the two confirmed 
cases in the Jordanian cluster died, we estimated that the 
expected cluster size could be 2/0-38, which is roughly 5. 
Hence, we considered three additionally probable cases 
next to the two confirmed ones in the Jordanian cluster. 
We also did a sensitivity analysis of our results, varying 
the size of the Jordanian cluster between two and 13 
(appendix). 

The partition of clusters into transmission trees is 
often ambiguous when the data are insufficient. We 
dealt with this limitation in several ways. First, we used 
published epidemiological investigations of MERS-CoV 
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clusters (ie, the Jordanian,” UK,* French,” Al-Hasa,” 
and Riyadh* clusters) to provide information about their 
tree structure. Second, we regarded the MERS-CoV 
clusters detected outside the Middle East as trees 
generated by the patient travelling outside the Middle 
East. Third, we made assumptions on the occurrence of 
transmission events on the basis of epidemiological 
links between cases to investigate the clusters that 
remained otherwise unresolved. We regarded cases as 
epidemiologically linked if the patients were either close 
family members or had documented contact and, 
furthermore, if the dates of onset were separated by the 
minimum incubation period (ie, 2 days”). Hence, we 
did not consider infectiousness outside the symptomatic 
period, since no data support it at present, and all 
transmission data available so far are compatible with 
transmission during the symptomatic phase.* This 
feature was also true for SARS-CoV? 

We investigated two different scenarios for partition- 
ing the remaining clusters. Our first scenario assumes a 
large number of index patients per cluster, each generat- 
ing a small tree, which results in a high introduction 
rate of the virus into the human population and a 
moderate transmissibility. The two cluster partitioning 
rules of the first scenario are that a report of an 
epidemiological link between two infected individuals 
certifies that transmission occurred between these 
individuals, and that all unlinked individuals are 
sporadic cases. However, we note that two epidemio- 
logically linked cases might have become infected after 
common MERS-CoV exposure. Hence, the first rule is a 
modelling choice that can be refined as further data 
become available. According to scenario 1, we regarded 
two Saudi Arabian clusters of two and four cases (the 
cluster of four cases located in Riyadh”) as single trees 
since they consisted of family members. The cluster of 
five reported in Dammam was partitioned into three 
sporadic cases and one tree of two cases, since two of 
the patients shared a hospital room. We split the 
Jordanian cluster into three sporadic cases and a tree of 
two cases (patient 3 and his brother, patient 8).” 

Our second scenario assumes a small number of large 
trees, which results in a low introduction rate and raised 
transmissibility of the virus. In this case, we grouped the 
individuals in one cluster to obtain the maximum tree 
sizes compatible with the epidemiological investigations. 
Hence, in the second scenario the Dammam and 
Jordanian clusters were regarded as transmission trees 
each containing five cases. The distribution of tree sizes 
obtained for both scenarios is given in the table. 


Data analysis 

We analysed MERS-CoV transmission using the theory of 
homogeneous branching processes. The key element of 
this theory is the distribution of the number of cases 
caused by an infected individual. Its average, the basic 
reproduction number of the pathogen infection R,, has a 


fundamental role in the transmission dynamics. If R, is 
less than 1, then all transmission trees terminate— 
otherwise transmission trees might be infinite and the 
disease becomes an epidemic. Assuming a Poisson 
distribution of the number of secondary cases, we 
inferred the R, of MERS-CoV using Bayesian analysis 
(appendix). We also did two sensitivity analyses of our 
R, results. First, we varied the size of the Jordanian cluster 
between two (the number of confirmed cases) and 13 (the 
total number of confirmed and probable cases). Second, 
we investigated the consequences of the outstanding 
event of six secondary cases caused by patient C of the 
Al-Hasa cluster in the dialysis ward (appendix). 

We calculated the rate of MERS-CoV introductions into 
the human population from the estimated number of 
index patients during the period of data collection (table). 
We assigned corresponding CIs on the assumption that 
introduction events follow Poisson statistics (appendix). 


Role of the funding source 

The sponsors of the study had no role in study design, 
data collection, data analysis, data interpretation, or 
writing of the report. The corresponding author had full 
access to all the data in the study and had final 
responsibility for the decision to submit for publication. 


Results 

Overall, patients with confirmed MERS-CoV were aged 
between 2 and 94 years (median 56 years, IQR 41-68 -5) 
and most were men (44 [74%] of 61 patients with 
sex reported). The case-fatality ratio was 59-4% 
(95% CI 46-71; 38 of 64 patients). Advanced age was a 
significant risk factor for death (case-fatality ratio 76% 
[55-91; 19 of 25] for age =60 years vs 38% [22-56, 
13 of 34] for age <60 years; p=0-008). 

We obtained an R, of 0-60 (95% CI 0-42-0-80) for 
scenario 1 and 0-69 (0-50-0-92) for scenario 2. The 
number of transmission trees for each scenario provided 
the number of MERS-CoV introductions into the human 
population during the period of data acquisition. 
Hence, we obtained yearly introduction rates of 
22-3 (95% CI 13-0-31-5) and 17-1 (9-0-25-3) for 
scenario 1 and 2, respectively. We also obtained an R, of 


Scenario 1 Scenario 2 
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Data are number of transmission trees of each size (size shown in far-left column) 
in each scenario. 


Table: Transmission tree sizes resulting from our interpretation of 
Middle East respiratory syndrome coronavirus data 
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Figure 2: Mathematical modelling results 


(A) The average tree size versus R, as predicted by the theory of homogeneous branching processes. The dark blue and red dashed lines correspond to the values of R, 
that we calculated for scenarios 1 and 2, respectively. (B) Contour plot of the expected Middle East respiratory syndrome coronavirus (MERS-CoV) yearly incidence 
versus introduction rate into the human population and R,. The very dark blue region corresponds to yearly incidence estimates below 10. The dark red region 
corresponds to yearly incidence estimates above 320. All other solid colour regions correspond to yearly incidence estimates bounded by the values shown on the 
contours. The black square and red circle show the parameter sets of scenarios 1 and 2, respectively; the error bars represent the corresponding 95% Cls. As a visual 
aid, we have shaded the region comprising the parameter sets compatible with scenarios 1, 2, and the in-between area compatible with intermediate scenarios. 

(C) The probability that R, exceeds 1 versus the size of the next MERS-CoV transmission tree; the horizontal dashed line corresponds to the 5% probability. (D) The 
probability that R, exceeds 1 versus the size of the next count of secondary cases of an index patient; the horizontal dashed line corresponds to the 5% probability. 


R,=basic reproduction number (the number of secondary cases per index case in a fully susceptible population). 


0-80 (0-54~-1-13) for prepandemic SARS. Repeating the 
analyses assuming the size of the Jordanian cluster to be 
between two and 13 or excluding the outstanding event 
of six secondary cases caused by patient C of the Al-Hasa 
cluster in the dialysis ward* showed the robustness of 
our results (appendix). As shown in figure 2A, the 
average tree size is very sensitive to the value of R,, 
particularly when R, is close to 1. 

We further studied the expected yearly incidence of 
MERS-CoV infections, calculated by multiplying the 
introduction rate into the human population by the 
average tree size, 1/(1—R,)." Figure 2B shows a contour 
map of the yearly incidence versus the introduction rate 
and R,. Non-linear dependence qualifies R, as the key 
parameter, particularly for values close to 1. On the same 
figure 2B, we plotted the parameters corresponding to 
our scenarios 1 and 2 and showed the region in the 
parameter space most likely to include the true 


parameters of MERS-CoV. This plot suggests that the 
possibility of MERS-CoV having an R, above 1 or less 
than 0-4, or a yearly introduction rate above 35, is very 
small. To assess the potential effect of future MERS-CoV 
outbreaks on our R, estimation, we analysed the change 
in the probability that R, exceeds 1 versus the tree size 
that would be observed next (figure 2C). In scenario 2, a 
tree larger than 98 cases would imply that the probability 
of having a pandemic MERS-CoV strain exceeds 5%. In 
scenario 1, reaching the 5% threshold would need the 
next observed tree to consist of 200 cases. 

When active contact tracing is operational, public 
health authorities might immediately stop viral trans- 
mission after the index patient has been diagnosed, and 
thus the available data would be just the number of 
secondary cases of the index patient before their isolation. 
This situation is addressed in figure 2D in which, in 
scenario 2, the 5% probability threshold is reached if the 
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next index patient caused eight secondary cases. In 
scenario 1, the 5% threshold is reached if the next index 
patient caused 14 secondary cases. 


Discussion 

Our analysis suggests that MERS-CoV has not reached 
epidemic potential—ie, R, was less than 1. Although R, 
estimates of prepandemic SARS-CoV and MERS-CoV 
under our second scenario are close, we cannot conclude 
that MERS-CoV will follow a similar path toward 
pandemic spread. Despite phylogenetic similarities, the 
two viruses have distinct biology, such as the use of 
different human receptors.°” SARS-CoV adaptation to 
human beings took just several months, whereas MERS- 
CoV has already been circulating for more than a year in 
human populations without mutating into a pandemic 
form. However, the speed of adaptation depends on 
several variables, some of which might widely differ 
between China and the Middle East (eg, animal and 
human density, baseline preventive measures, and, most 
importantly, the prepandemic R,)." Because of the very 
small number of MERS-CoV sequences available, non- 
synonymous mutations suggesting viral adaptation could 
not be yet established. Hence, comparison between SARS- 
CoV and MERS-CoV should remain cautious. 

Figure 2B displays the yearly number of cases as a 
function of two variables, the yearly introduction rate and 
R,. Control of this new zoonosis could therefore be 
exercised by addressing these two components. 

First, control could target a reduction in the rate of 
MERS-CoV introductions into the human population. 
However, the animal host responsible for MERS-CoV 
introductions is not yet known. This search should take 
high priority, focus on Middle Eastern countries where 
most cases occurred, and address the broad variety of 
species known to host coronaviruses, including bats, 
birds, mice, dogs, pigs, and cattle.” Reports have already 
pointed to sick goats and camels in the vicinity of patients 
as a potential source of MERS-CoV.” 

Second, control could focus on a reduction in the infec- 
tious period of infected individuals through improved 
surveillance, rapid diagnosis, and isolation. When done 
effectively, active contact tracing could reduce the 
number of cases to the number of MERS-CoV intro- 
ductions (ie, index patients) and the secondary cases they 
might cause before isolation.” We note that R, might 
also increase with the population density and is affected 
by the community age and contact structure. Large 
gathering events or travel of patients from the Middle 
East to densely populated areas might result in outbreaks 
with increased R,. 

Our analysis depends on the quality of the surveillance 
systems.” In view of the tendency for identification of large 
clusters at the beginning of an epidemic, and insufficient 
data for partitioning clusters into smaller transmission 
trees, our tree-size data analysis might overestimate R,. 
Hence, this bias would not affect our conclusion that 


MERS-CoV does not have pandemic potential and is likely 
to be minimised now that public health surveillance 
systems are on alert. Also, timely control measures after 
diagnosis of a new patient will trim the natural pattern of 
transmission. Still, the virus will always have a short 
window of opportunity for transmission before the iso- 
lation of each new index patient. For this reason, we 
considered the situation in which the next available 
datapoint would be the count of cases secondary to the 
index patient, and provided a method that might be useful 
in places with effective surveillance. 

Another important issue is the possibility of asympto- 
matic and mild infections. A confirmed case with mild 
symptoms, for which admission to hospital was not 
needed, occurred in the UK.‘ If asymptomatic and mild 
infections constituted an important fraction of the total, 
R, might be higher than our estimate. In the future, 
serological analysis of contacts will allow a fair estimation 
of asymptomatic infections. Serological tests for MERS- 
CoV are under development, but have not yet been 
deployed. Retrospective analysis of sera of contacts of 
SARS-CoV patients showed a very small (0-2%) pro- 
portion of asymptomatic infections.* Of note, serum 
samples of 2400 control patients at the Dr Soliman 
Fakeeh Hospital in Jeddah, Saudi Arabia, where the first 
patient was diagnosed, did not detect MERS-CoV in an 
immunofluorescence assay that was strongly positive 
with the serum of the patient, suggesting that MERS- 
CoV was not circulating undetected in the general 
population during the previous 2 years.’ 

The main issue for now is how to interpret the data on 
the size of future outbreaks that surveillance systems will 
detect. Most certainly, future data will add to our 
knowledge about viral transmissibility and narrow the 
Cls around the R, estimate. Our analysis suggests that if 
effective contact tracing detected eight or more secon- 
dary cases caused by the next index patient, then, 


Panel: Research in context 


Systematic review 

A PubMed search on June 23, 2013, with the terms 
“HCoV-EMC", "MERS-CoV”, and “novel human coronavirus” 
identified 24, 14, and 46 reports, respectively, linked to 
Middle East respiratory syndrome coronavirus (MERS-CoV) 
infection. 18 reports contained detailed epidemiological 
information about MERS-CoV cases, starting with the Sept 27, 
2012, description of two cases by Danielsson and colleagues.” 
Articles relevant to our paper are cited in the text. 


Interpretation 

Our report provides the first estimation of the basic 
reproduction number R, for MERS-CoV, with evidence that 
the virus has not yet reached pandemic potential. Enhanced 
surveillance, active contact tracing, and vigorous searches for 
the MERS-CoV animal host and transmission route to human 
beings should be urgently prioritised. 
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under scenario 2, this would suggest, together with all 
knowledge gathered so far, that there is a more than 5% 
chance that the MERS-CoV R, is larger than 1—ie, that 
MERS-CoV might have pandemic potential. This finding 
seems unlikely because the largest count of secondary 
cases attributed to a patient is seven (six in a dialysis 
ward and one outside) and is an outlier. 

In conclusion, our analysis confirms the importance of 
enhanced surveillance of MERS-CoV infection and an 
active search for its animal host, particularly in the 
Middle East where MERS-CoV infection is most prevalent 
(panel). Close monitoring of cases and contact tracing 
are of high priority to limit transmission, gather high- 
quality data to update the R, estimates, and decrease 
opportunities for viral adaptation to interhuman 
transmission. One of the main lessons of the SARS 
pandemic has been that early control of the virus (while it 
was still confined to southeast China) might have 
prevented its global spread. 
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